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The dynamics of escape from an attractive state due to random perturbations is of central interest 
to many areas in science. Previous studies of escape in chaotic systems have rather focused on the 
case of unbounded noise, usually assumed to have Gaussian distribution. In this paper, we address 
the problem of escape induced by bounded noise. We show that the dynamics of escape from an 
attractor's basin is equivalent to that of a closed system with an appropriately chosen "hole" . Using 
this equivalence, we show that there is a minimum noise amplitude above which escape takes place, 
and we derive analytical expressions for the scaling of the escape rate with noise amplitude near the 
escape transition. We verify our analytical predictions through numerical simulations of two well 
known 2-dimensional maps with noise. 



The escape of trajectories from attracting sets due to 
the effect of noise has been a central issue in various 
branches of science for a long time. From a fundamen- 
tal perspective, the study of the dynamics of escape in- 
cludes the fundamental work by Arrhenius on chemical 
reactions [1] , passing through ideas of Kramers from the 
forties of last century [21, the escaping on chaotic dy- 
namics [SHSj to very recent work on Statistical Mechan- 
ics [7]. Notwithstanding this long history, an important 
case has mostly been neglected, namely the dynamics 
under bounded noise. In fact, the vast majority of the- 
oretical works in this area is heavily dependent on the 
assumptions of unbounded noise, almost always assumed 
to have a Gaussian distribution [31 [71 |H] ■ Thus, they are 
not applicable to other cases. In particular, very little 
is known about the dynamics of systems with escape in 
the presence of bounded noise. As for many applica- 
tions a bounded perturbation is arguably more realistic 
than unbounded ones, this is an important gap in our 
understanding. For example, in Neuroscience one may 
be interested in the minimum energy for bursting to take 
place [5] , in Geophysics one may consider the overcoming 
of some potential barrier just before an earthquake IlOj, 
or the critical outbreak magnitude for the spread of epi- 
demics [11], among others. 

In this paper, we investigate escape in dynamics per- 
turbed by bounded noise using a new approach. We de- 
scribe the noisy dynamics in terms of a family of random 
maps, whose iteration gives rise to a discrete Markov pro- 
cess with transition probability supported in the neigh- 
bourhood of the points generated by the iteration of the 
deterministic dynamics. The dynamics near an attractor 
is then described using the formalism of conditionally in- 
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variant measures. We show that the dynamics of the 
escape is determined by the measure of a subset Iq of 
the phase space, which we call conditional boundary. It 
consists of those points lying close enough to the basin 
boundary such that they are subject to escape under the 
effect of random perturbation in one iteration of the map, 
and intersecting the support of the conditionally invari- 
ant measure of the system. Therefore, the dynamics is 
mapped onto the dynamics of a closed system with a hole 
— Iq being the hole. We show that there is a minimum 
noise amplitude for the escape to take place, which is 
the critical noise amplitude that makes Iq non-empty. 
We show that the mean escape time is determined by 
the measure of the conditional boundary, and use this to 
derive a power law relation between the average escape 
time (r) and the noise amplitude i^, for ^ close to (and 
higher than) 

(1) 

where a depends on the dimension of the system. We 
show that for dimension two, a = 3/2, and we verify 
that this prediction is correct by comparing with the re- 
sults of numerical simulations for two maps from different 
families. This result is independent of any particular sys- 
tem, and holds universally for bounded noise. Equation 
^ is in contrast with the case of Gaussian noise, where 
an exponential scaling is observed [THE]- Furthermore, 
although the noise-induced escape from attractors may 
seem to be very different from that of systems undergo- 
ing a boundary bifurcation, we show that our approach 
allows us to establish a connection between these two pro- 
cesses, and to explain why we find the same time scaling 
in both cases. 

To get started, first consider a deterministic dynam- 
ics Xn+i — .f{xn) given by the iteration of the map /, 
a smooth function whose inverse is differentiable, in the 
phase space of the system M (i.e., the iteration of a dif- 
feomorphism f : M ^ M). Our main focus will be on 
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invariant subsets A of M which attract their neighbour- 
ing points, that is, /"(x) tends to A as n — cx); these 
are the attractors of the system. The basin of attraction 
of A is the open set W''{A), the set of points eventually 
coming close and converging to A. The next ingredient is 
the boundary of the basin of attraction, a zero Lebesgue 
measure ergodic component of the phase space which we 
denote by d. 

We shall consider a random perturbation of the deter- 
ministic system introduced above [121 US], so that our 
perturbed system is described by a family of random 
maps, that in our context can be written as 

F{x,)^f{x,)+s,, (2) 

with \\sj\\ < ^, where Sj is the vector of random noise 
added to the deterministic dynamics at the iteration j, 
and ^ is its maximum amplitude. In this way, is a con- 
tinuous application [25]. We illustrate this in Fig. [TJa). 
As it is shown in the picture, the perturbed dynamics 
can be thought as follow. Image we iterate the point Xj 
by the deterministic system /. Then, let say that at the 
moment that we take the f{xj) to evolve our dynamic 
again, we make a small error given to our limited preci- 
sion. Nevertheless, we can assure that the error is always 
less than ^. Then we ask whether the attractors for the 
system with no error are still attractors when a small er- 
ror is considered. If so, how large can our error be such as 
the attractors will still be preserved? Above this thresh- 
old, how does the escape of orbits scale? These natural 
extensions of deterministic processes are exactly the sort 
of problems we shall be interested in. 

From the probabilistic perspective, the idea of orbits 
converging to some attractor is represented by the con- 
cept of physical (or SRB) measures [131 [T3] . Suppose ini- 
tially that we compute the time average, ^ j=o^ (^) ' 
along the orbit for a given initial condition as n evolves. 
This quantity, the time average, is expected to converge 
to some invariant value as the orbit approaches the at- 
tractor. If we extend it for a large number of initial condi- 
tions, then, the time averages that we compute along dif- 
ferent orbits are also expected to converge to an invariant 
value as these orbits approach the attractor. It turns out 
that such value defines the so-called SRB measure, which 
characterises the attractor. The set of initial conditions 
that we chose for computing such time averages along the 
orbits are exactly what we call the basin of the measure, 
that we represent by Since we want to statistically 

characterise the behaviour of the attractor, it is desirable 
the set of initial conditions for which we compute such 
time averages to have positive Lebesgue measure. That 
is, there is a non negligible number of initial conditions 
whose time averages along the orbits converge to such 
invariant quantity and thus characterise the attractor. 
Therefore, we say that a measure is physical if its basin 
has positive Lebesgue measure. More precisely, the basin 
^{fi) of such measure is the set of points whose time av- 
erages along the orbits weakly converge to the space av- 
erage, Hm„^oo J2]Zo V^if^i^)) = / Vdf^, where (p rep- 
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FIG. 1: (a) Illustrative picture of the perturbed dynamics 
with amplitude of noise ||ej|| < (b) A basin of attraction 
W''{A) and its basin boundary d (dashed hue). We illustrate 
that the iteration z i— >■ f{z), from the point z initially in 
W (A) brings the orbit within a distance ^ from the boundary. 
Therefore, the random perturbation applied to f{z) with some 
|£|| < ^ could push the random orbit outside the basin. On 
the other hand, for the point x G W''{A) the iteration x i— >■ 
f(x) brings it farther than the maximum perturbation ^ away 
from the boundary d. Therefore, in our illustration z £ Ig 
but X ^ Ig. 

resents our measurements of an observable |26j . Here, 
a physical measure of a set A C M is an invariant er- 
godic probability measure /i supported on A. Therefore, 
such measures are closely related to the so-called natural 
measures. Note that as a consequence of the Birkhoff's 
ergodic theorem, 

^Ji{^^{^l)) = 1. (3) 

In words, Eq. [3] just says that we are normalising such 
measure on the total set of initial conditions whose time 
averages converge to the invariant one. For most known 
cases the basin of the physical measures supported on the 
attractors coincides with the basins of the attractors [13] , 
so we assume here that, W^^(A) — ^(/i). This is known 
as the basin property. 

Despite the complications introduced by the presence 
of noise in the dynamics, for continuous random maps 
with small amplitude of perturbation an invariant prob- 
ability measure is guaranteed to exist. This is due to 
the Krylov - Bogolubov Theorem |15j , which ensures that 
every continuous application in a compact measurable 
space has an invariant probability measure. For ran- 
domly perturbed systems, we expect that the perturbed 
orbit will densely fill up the neighbourhood of the attrac- 
tor. This property is called random transitivity (12j . and 
under general conditions it can be shown that these mea- 
sures are ergodic |12j . As an example, we can think of an 
attractive fixed point for a deterministic system, and the 
perturbed version of that system. In the perturbed dy- 
namics, there is a density of probability around the fixed 
point. When we increase the amplitude of the noise, the 
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density of probability becomes more spread around the 
original fixed point. Conversely, decreasing the ampli- 
tude of noise, the support of the invariant measure tends 
to the point attractor. 

Now consider a subset Iq of the phase space neigh- 
bouring the basin boundary, defined by the set of points 
X whose f{x) is within a distance ^ from some point in 
the boundary 9, 



(a) 
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where B{f{x),^) is the ball of radio ^ around f(x). This 
is illustrated in Fig. [ijb). In addition to being close 
enough to the boundary, another important condition for 
the escaping process to take place is that the intersection 
of the support of the invariant physical measure of the 
system under stochastic perturbation, supp/zj, overlaps 
the set Id- That is, the density of probability around the 
attractor needs to be spread over a region close enough 
to the boundary. We thus define the set Iq by 



Iq^ IqC] supp/i^. 



(5) 



Because the noise is bounded, for very small noise am- 
plitudes Iq = ^. But as the amplitude of the noise is 
increased, we expect that for a certain critical amplitude 
^ = ^c, we have 7^ 0, and escape takes place for any 
^ > ^c- We call Iq the conditional boundary. The im- 
portance of Id{C) stems from the fact that it represents 
the set of points which one iteration of the map / can po- 
tentially send close enough to the boundary d, such that 
a random perturbation with amplitude ^ may send them 
out of the basin of attraction. The dynamics of escape is 
thus governed by this set, and it can be understood as a 
"hole" which sucks orbits from the basin if they land on 
it. 

Due to property ([s]), we have that /g ^ for ^ > 
and thus > 0. We can therefore think of the system 

with ^ > as a closed system with a hole (or leak) [T^] , 
where Iq plays the role of the hole. The idea of a system 
with a hole has been used in different contexts before, 
for example |16j . Furthermore, note that the measure 
is not in fact invariant because of the loss caused 
by escape. However, it is possible to describe such escape 
problem in terms of conditionally invariant measure [161 
117] . We say that a probability measure /ic is conditionally 
invariant with respect to F if 
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for every measurable subset X C A, where A C M is 
a non-invariant region of the phase space M. The con- 
ditional measure is defined in a way that, for each it- 
eration, when the set A loses a fraction of its orbits to 
the hole, we renormalise its measure by what remains in 
A. It is defined in terms of pre-images (or more gener- 
ally using pushforward measures [16].) In our context, 
the hole is the set Iq, and therefore fic{Id) > is pre- 
served by the dynamics because of the compensation fac- 
tor fj,c{F~^{A)) [roi[TH]. Due to the random transitivity 
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FIG. 2: We illustrate the basin boundary d, the region Ig 
whose edge is the dashed line /3, and the boundary of the 
support of the invariant measure for the perturbed system 
(the parabolic line). In (a), ^ < ^c, thus Ig n supp/xj = 0. In 
(b), (, — the supp/ij is tangent to the limit of Ig. In (c), 
£, > (,c- As a consequence, Ig n supp therefore, Ig ^ ^ 
(shadowed area), what implies fi{Ig) > 0. 



of the conditionally invariant domain, we expect no par- 
ticular dependence on the density of points in A. In this 
case, a random trajectory diffuses through the support 
of /ic, until it eventually comes inside Iq. Once in this 
set, there is a probability that it will permanently escape. 
Indeed it is the measure of Iq which controls the escape 
rate. Because we are interested in the regime of ^ ^ ^c, 
hence small leaks, we can assume fT51. 



Hc{Id) = l^ih) > 0. 



(7) 



Since there is a certain probability that a particle escapes 
if it falls into Iq, we expect from Kac's Lemma [19] that 
the average escape time satisfies 
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1 



(8) 



Rigourously proving the existence of and calculating 
fJ.{lQ) is not an easy task. We use here a heuristic ap- 
proach to obtain the scaling of ^(/g) for ^ close to [20] . 
For simplicity we focus on the 2-dimensional case. For 
C < Cc; we expect the probability distribution to be rela- 
tively concentrated on the attractor, thus having support 
with Lebesgue measure (area) smaller than that of the 
basin of attraction. As a first approximation, we image 
the edge of the support as being a smooth closed curve. 
As ^ is increased, the mean radio of the distribution grows 
but it needs to be less than that of the basin of attraction, 
otherwise we would have the invariant measure supported 
outside the basin of attraction. When ^ ^ ^c, we picture 
it as a generic intersection of two curves of different ra- 
dios. Therefore, the portion of the edge of the support 
of the invariant measure lying within Ig is locally well 
approximated by a parabola (curved line in Fig. [2]) . For 
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C ~ Cc, we have Ig ^ 0, thus fidld) > 0. Recalhng Eq.j?) 
it gives us that > 0. 

To estimate the measure of the hole, define 



(9) 



The top curve encompassing the shadowed area in 
Fig. [2jc) is then weU approximated by 



(10) 



in the appropriate units. S intersects /3 in two points, 
namely and y/fj. Because we assumed the basin 

property to hold, in other words, we normalised the mea- 
sure and assumed the basin of the measure to be equal to 
the basin of attraction of the deterministic system, Eq. [3] 
also tells us that the probability measure of the hole is 
proportional to its Lebesgue measure, the area, encom- 
passed for /3 and S within the interval [— y^, y/fj] — the 
shadowed area in Fig. [2]jc). We can calculate it as 



, and thus 
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M^a)«(e-ec)^/'. 



Applying the Kac's Lemma, we obtain Eq. [T] with a = 
3/2. Note that although we use the approximation to de- 
scribe the boundary d locally as a smooth curve, it might 
actually be fractal. Therefore, if the random orbit falls 
into Ig, there is only a probability that it will escape due 
to the fractal property of d. This fact is subtly incor- 
porated by Eq. [8] in the proportionality rather than the 
equality to the inverse of the mean escape time. 

In order to check this prediction we numerically ob- 
tained the scaling of the distribution of escaping times 
with amplitude of noise for two distinct 2-dimensional 
systems. The first perturbed systems we have chosen was 
the randomly perturbed single rotor map [U], defined 
by 
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J I yj{mod2'K) 
u)yj + 4sin(a;j 4 



(12) 

where x € [0, 27r], and y £ M, and represents the dissi- 
pation parameter. As a second testing system, we have 
chosen the perturbed dissipative Henon map, in the form 



G 



_ ( 1.06x2 - (1 - y)y.^ 



V] 



(13) 



where, x and y are real numbers and again, represents 
the dissipation parameter. We used v = 0.02 for both 
maps, as for such value they present very rich dynam- 
ics [22^. We also assumed, for the purposes of this nu- 
merical experiment, the noise to be uniformly distributed 
in each variable; but we stress that this is just a numeric 
convenience. For each map, we computed the time that 
random orbits took to escape from their respective main 
attractors for a range of noise amplitudes. In each case. 
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FIG. 3: (Colour online) The inverse of the mean escape time 
scaling with amplitude of noise for the Map ( 12 1 - black circles 
- and for the Map ( 13 1 - blue squares. For each map, the 
values of mean escape time were obtained by iterating 10^ 
random orbits for each value of ^. The dashed lines show the 
expected scaling (^ — ^c)^''^ and the thick continuous lines 
show the best fitting for the Rotor map, a ~ 1.7, and for the 
Henon map, a « 1.6. In the insets, we show the transition. 
For ^ < ^c, we have (T) = oo, the random orbits do not 
escape, therefore 1/{T) = 0. For ^ > ^c, the escape time 
scales as Eq. |[l]). For the used parameters ~ 0.086 ±0.006 
for the Rotor map - inset (a) - and = 0.021 ± 0.002 for the 
Henon map - inset (b). 



the mean escape time was obtained for 10'^ random orbits 
for each value of ^. The results are shown in Fig. [3] For 
the parameter used here, we obtained = 0.086 ± 0.006 
for the perturbed Rotor map and — 0.021 ± 0.002 for 
the perturbed Henon map, what is shown in the insets. 
In both cases, we obtained a good agreement between our 
simulations and the predictions of our theory for a range 
of decades. An important remark is regarding the preci- 
sion of the ^c- For the one dimensional case, for example, 
it has been proved that similar power laws in the unfold- 
ing parameters are in fact lower bounds for the average 
escape time scale |17j . Therefore, even from the numer- 
ical perspective, it is difRcult to accurately estimate the 
value of ^c- Indeed, in our case we observe when increas- 
ing ^ near some transient irregular bursts regime before 
the escape phenomena becomes robust. For a number of 
initial conditions, we have thus a distribution of values 
of critical noise around ^c, that is expected to become 
sharper as the number of initial conditions is increased. 
As a direct implication, the exponent obtained in our nu- 
merical simulations also varies within some range. For ex- 
ample, for the Rotor map, if we choose = 0.086-1-0.006, 
we obtain a = 1.5 and for — 0.086 — 0.006, we have 
a — 2.0. For the Henon map, we obtained a = 1.3 for 
£_c = 0.023 and a = 2.1 for £_c = 0.019. 

Note also that, in principle, a much larger number of 
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random orbits would be necessary for one to be able to 
observe a "perfect" power law. This is because most of 
the theoretical arguments used here, such as the conver- 
gence of time averages, are obtained in the asymptotic 
limit. Furthermore, we notice that for "large" values of 

— ^c), meaning large holes, our simulated results differ 
appreciably from our theoretical prediction. This is due 
to the fact that /ic(^a) = is valid only for small 

leaks [TH]. In addition, for large amplitude of noise, 
the dynamics is totally dominated by the noise, and is 
not well described as a small perturbation around the 
deterministic motion. We note here that the exponen- 
tial distribution reported in ^ for a particular case of 
bounded noise was obtained for values of ^ much greater 
than the critical amplitude which is an outside the 
range of validity of our theory. 

As last consideration, we want to call attention to the 
case where noise is applied on bifurcating systems [71 1231 
124] . In our approach we consider rj to be increasing. In 
Fig. [2] this corresponds to moving the parabola upwards 
until it intersects the /3 curve, when the transition to 
escape takes place. We can easily see that we should 



expect an equivalent transition to escape by moving /3 
instead, as a result of changing a bifurcating parameter 
whilst keeping the noise amplitude constant. Therefore, 
the exponent obtained for the case of bifurcating systems 
is expected to be the same as ours, and this is indeed the 
case [a EH Ell- 
in conclusion, we have shown that the problem of es- 
caping orbits from attractors due to the effect of bounded 
random noise can be thought as a closed system with a 
hole. We identify the subset of the phase space that is 
responsible for the escape and acts as a hole, and show 
that the measure of this set determines the escape rate. 
We have shown that there is a critical amplitude of noise 
in order to such escape happens. When the amplitude 
is just above this critical value, we derived a universal 
power-law relation of escape time with respect to the 
amplitude of noise, in contrast with the case of Gaussian 
noise. 

The authors a grateful to the anonymous referee for 
valuable suggestions to increase the value of this paper. 
AM and CG have been supported by the BBSRC, under 
grants BB-F00513X and BB-G010722. 



[1] P. Hanggi J. Stat. Phys. 42, 105 (1986). 
[2] H. A. Kramers Physica (Utrecht) 7, 284 (1940). 
[3] P. Grasberger J. Phys. A 22, 3283 (1989). 
[4] S. Kraut, U. Feudel, and C. Grebogi, Phys. Rev. E 59, 
5253 (1999). 

[5] S. Kraut, and U. Feudel, Phys. Rev. E 66, 015207 (2002). 
[6] S. Kraut, and C. Grebogi, Phys. Rev. Lett. 92, 234101 
(2004). 

[7] J. Demaeyer, and P. Gaspard Phys. Rev. E 80, 031147 
(2009). 

[8] P. D. Beale Phys. Rev. A 40, 3998 (1989). 
[9] N. Nagao, H. Nishimura, and N. Matsui, Neural Process- 
ing Lett 12, 267 (2000); S. J. Schiff and K. Jerger and D. 
H. Duong and et al., Nature 370, 615 (1994). 

[10] O. Peters, and K. Christensen, Phys. Rev. E 66, 036120 
(2002); P. Bak, K. Christensen, L. Danon, and T. Scan- 
ion, Phys. Rev. Lett 88, 178501-1 (2002); M. Anghel, 
Chaos SoUt & Frac 19, 399 (2004). 

[11] L. Billings, E. M. BoUt, and I. B. Schwartz, Phys. Rev. 
Lett 88, 234101 (2002); L. Billings, and I. B. Schwartz, 
Chaos 18, 023122 (2008). 

[12] V. Araujo, Ann. Inst. Henri Poincare, Analyse non 
lineaire 17, 307 (2000). 

[13] C. Bonatti, L. J. Diaz, and M. Viana, Dynamics beyond 
uniform hyperbolicity Springer, Berlin (2005). 

[14] J. -P. Eckmann, and D. Ruelle, Rev. of Modern Phys. 57, 
617 (1985). 

[15] R. Mane, Ergodic theory of differentiable dynamics 

Springer Verlag (1987); 
[16] G. Pianigiani, and J. A. Yorke, Trans. Am. Math. Soc. 

252, 351 (1979); M. F. Demers, and L.-S. Young, Non- 



lineanty 19, 377 (2006); A. E. Motter, and P. S. Letelier, 
Phys. Lett. A 285, 127 (2001); M. A. Sanjuan, T. Horita, 
and K. Aihara, Chaos 13, 17 (2003); L. A. Bunimovich, 
and C. P. Dettmann, Phys. Rev. Lett. 94, 100201 (2005); 
Europhys. Lett. 80, 40001 (2007); 
[17] H. Zmarrou, and A. J. Homburg, Ergod. Th. & Dynam. 
Sys. 27, 1651 (2007); Discrete Cont. Dyn. Sys. BIO, 719 
(2008). 

[18] E. G. Altmann, and T. Tel, Phys. Rev. Lett. 100, 174101 
(2008); Phys. Rev. E 79, 016204 (2009). 

[19] M. Kac, in Probability and Related Topics in Physi- 
cal Sciences, Intersciences Publishers, New York (1959), 
Chap IV. 

[20] Z. Liu, Y-C Lai, L. Billings, and I. Schwartz Phys. Rev. 

Lett. 88, 124101 (2002). 
[21] G. M. Zaslavskii, Phys. Lett. A 69, 145 (1978); B. 

Chirikov, Phys. Rep. A 52, 265 (1979). 
[22] C. S. Rodrigues, A. P. S. de Moura, and C. Grebogi, 

Phys. Rev. E 80, 026205 (2009). 
[23] E. Ott, Chaos in Dynamical Systems 2ed., Cambridge 

(2002). 

[24] J. C. Sommerer, E. Ott, and C. Grebogi, Phys. Rev. A 
43, 1754 (1991). 

[25] The iteration of random maps can also be seen as an 
associated discrete time Markov process having a family 
{p^{-\x) : X £ U,^ > 0} of transition probabilities, where 
every p^{-\x) is supported on U, a ^-neighbourhood of 
/(^)- 

[26] More generally, for every continuous function (/p : A/ — >■ R. 



